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The Harper-Hofstadter model provides a fractal spectrum containing topological bands of any 
integer Chern number, C. We study the many-body physics that is realized by interacting particles 
occupying Harper-Hofstadter bands with \C\ > 1. We formulate the predictions of Chern-Simons 
or composite fermion theory in terms of the filling factor, v, defined as the ratio of particle density 
to the number of single-particle states per unit area. We show that this theory predicts a series of 
fractional quantum Hall states with filling factors v = r/(r|C|+ 1) for bosons, or v = r/(2r|C|-|-l) for 
fermions. This series includes a bosonic integer quantum Hall state (bIQHE) in \C\ = 2 bands. We 
construct specific cases where a single band of the Harper-Hofstadter model is occupied. For these 
cases, we provide numerical evidence that several states in this series are realized as incompressible 
quantum liquids for bosons with contact interactions. 


Recently, there has been much progress towards exper¬ 
imental realizations of topological flat bands, such as by 
light-matter coupling in cold gases [1-7] or via spin-orbit 
coupling in condensed matter systems [8-10]. These sys¬ 
tems provide novel avenues for exploring fractional quan¬ 
tum Hall physics in new settings where lattice effects play 
important roles [8-19]. Furthermore, these “fractional 
Chern insulators” generalize the fractional quantum Hall 
states of interacting particles in continuum Landau levels 
to lattice-based systems. 

In cases where the underlying topological band has 
unit Chern number, C = 1, the states can be continu¬ 
ously connected to conventional fractional quantum Hall 
states in the continuum Landau level [20]. However, if 
the band has Chern number, C , of magnitude greater 
than 1, no such continuity is possible. The fractional 
quantum Hall states have features that are particular 
to the lattice structure. The appearance of fractional 
quantum Hall states for bands with |C7| > 1 has been 
demonstrated in various lattice models, with unit cells 
that contain multiple states in the form of distinct sub¬ 
lattices, or in terms of internal degrees of freedom such 
as spin or color [see, e.g. 21]. In particular, this has led 
to the proposal of states at filling factors v = 1/(|C| +1) 
for bosons [22-25]. 

In this Letter we show that this physics of interact¬ 
ing particles in the novel Chern bands can be captured 
within the Harper-Hofstadter model, in which a mag¬ 
netic unit cell arises naturally without additional inter¬ 
nal degrees of freedom. This model leads to a complex 
energy spectrum as a function of flux n q i = /H>o per 
plaquette. The low energy bands can have Chern num¬ 
bers larger than one. We show that they can realize the 
sequences of fractional Chern insulator states for \C\ > 1 
discussed for other models, providing an interpretation 
of these states in terms of the composite fermion con¬ 
struction on a lattice [11, 14], Based on these insights, 
we identify another sequence of fractional Chern insula¬ 
tor states with filling factors v = r/(r\C\ + 1). We show 
numerical evidence for this sequence from exact diago- 
nalization studies. 

The Harper-Hofstadter model has recently been re¬ 
alized, at least for weakly interacting particles, in ex¬ 
periments on ultracold gases [2-4]. Further realizations 


have also been obtained for a triangular moire lattice 
in graphene flakes deposited on boron nitride [26-28]. 
Our results demonstrate that, under suitable conditions 
(particle density, flux density and temperature), these 
systems have the possibility to explore a wide range of 
the novel physics of fractional Chern insulators. 

The rich physics of charged particles in a two- 
dimensional plane subjected to a perpendicular magnetic 
field B and a periodic potential, or physically equivalent 
models emulating this scenario, arises from its two com¬ 
peting length scales: the magnetic length £ 0 = ^Jh/eB 
and the lattice scale a. Hence, the problem brings to 
play the commensurability of these two scales, as first 
analyzed by Harper [29]. The resulting fractal structure 
of the single-particle spectrum was revealed by Azbel [30] 
and the full spectrum solved numerically and character¬ 
ized by Hofstadter [31]. As shown by Wannier [32], the 
Azbel-Hofstadter recursion relations imply that the total 
number of states per unit area n s below each gap varies 
linearly with the flux density n^, and further that their 
relationship is described by the Diophantine equation 


n s = Cn^ — D, C,DeZ. (1) 


The work by Streda [33], as well as Thouless et al. [34, 
35] explains the physical relevance of Wannier’s result. 
For noninteracting fermions, one expects incompressible 
states at density n = n s , with Hall conductivity given 
by Streda’s formula as 
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( 2 ) 


Thouless et al. derive the Hall conductivity by direct 
calculation from a Kubo formula, obtaining the integer 
quantization from the topological nature of the resulting 
expression, namely that C =^Ci where C; is the Chern 
number of the i-th occupied band [34], Thus, the integer 
C in Eq.(l) is seen to be the net Chern number of the 
bands contributing to the n s states below the energy gap. 
Solutions to (1) exist for any C in which n s arises from 
a single band (see below), establishing the presence of 
bands of any Chern number in the Hofstadter spectrum 
(albeit with rapidly decreasing gaps for large C ) [36, 37]. 

To realize the nontrivial Chern bands of the Harper- 
Hofstadter model, it is sufficient to create a tight-binding 




2 


lattice system with complex hopping elements between 
nearest neighbor sites, a fact exploited for realizations of 
the model in cold gases [2, 3, 6 ]. From the analysis of 
the magnetic translation group [38, 39], it follows that at 
a rational flux density n$ = p/q , the Harper-Hofstadter 
Hamiltonian admits a periodic representation on a mag¬ 
netic unit cell (MUC) comprising q sites, i.e., an area en¬ 
closing an integer number of flux quanta and an integer 
number of plaquettes of the lattice. The single-particle 
Hamiltonian then takes the tight-binding form 

ft s p = ~y] Uje^ a]ai + h.c., (3) 

i,3 


in which the phases <f>ij are invariant under translations 
of the MUC. (The choice of the MUC and the vector po¬ 
tential A make up the remaining space of gauge choices.) 
One can therefore consider the q sites of the MUC as 
sublattices a = 1 , of a general q -site tight-binding 
model and solve via Bloch’s theorem. Note that the ori¬ 
gin of the MUC can be chosen on any site of the lattice, 
so the problem has an additional g-fold symmetry. We 
give an analysis of the single-particle properties in the 
supplementary material [40], taking care to respect this 
symmetry [41]. 

To search for strongly-correlated phases, it is useful 
to identify situations in which there is a manifold of 
low-energy single-particle states (one band, or several 
closely spaced bands) that is well separated from higher- 
energy bands. For now, we focus on the case in which 
this manifold is a single band, with a large gap to the 
next band. The largest gap in the Harper-Hofstadter 
spectrum corresponds to the lowest Landau level, with 
C = 1, D = 0. Here, we seek more general states, and 
consider the next-largest gaps found at the first level of 
the Hofstadter hierarchy. These appear in the vicinity of 
cell boundaries close to the simple rational flux densities 
ric/, = 1/Q, (Q > 1), where the energy bands become 
exponentially flat in terms of both their energy disper¬ 
sion and their Berry curvature. Hence, these bands are 
well suited to support incompressible fractional quan¬ 
tum Hall states [13, 14, 42, 43]. The gaps near points 
(n s ,ri 0 ) = (0, Q^ 1 ) are described by the Diophantine 
equation ( 1 ) with C = sQ, and D = sgn(C) = s, with 
s = ±1 for the bands at n<f, ^ 1 /Q. We are particularly 
interested in cases where we find a single band below this 
gap, i.e., where n s (n^, = p/q ) = 1 /q, and thus q = Qp—s , 
with corresponding flux densities 


ricj, 
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Below, we take p > 2 to ensure that the band belongs 
to the subcell nearest n$ = 1/|C|. We will also consider 
higher band gaps at the same flux densities, which can 
be seen as fractal replicas of the r-th continuum Landau 
level, for which 


n s = r(C?X 0 - sgnC), r £ Z\{ 0}, (5) 

and where rC ^ 0 for the flux densities n ^ IQ- 
We now discuss the many-body physics of interacting 
particles in the Harper-Hofstadter model, described by 


the Hamiltonian 

ft = ftsp + \ V(ri - Tj) : flifij :, ( 6 ) 

i,3 


with site labels i, j, and : h : denoting normal ordering of 
the density operators. Let us first review the predictions 
of Chern-Simons theory [ 11 ], or equivalently the lattice 
composite fermion picture [14], and translate these re¬ 
sults into the language used for the analysis of Chern in¬ 
sulators. The basic premise of this approach is that the 
interaction includes a sufficiently strong short-range re¬ 
pulsion in order to favor “flux attachment,” which keeps 
the particles at a distance from each other, thus minimiz¬ 
ing interaction energy. The composite fermion Ansatz 
translates this idea into a trial wave function of the form 
^triai(ri,...,rjv) = T’lem^ , j({u}) x ^cf({u}), where 
both the Jastrow factor T.j and the composite fermion 
wave function Tcf vanish when the positions of two par¬ 
ticles coincide, and "Plem denotes the projection onto 
the relevant low-energy manifold of single-particle states. 
For the case of bosons (fermions), one needs to attach an 
odd (even) number k of flux quanta to the particles, so 
as to obtain an effective problem of weakly interacting 
composite fermions (CF) experiencing an effective flux 
density n*^ relating to the externally applied flux via 

ritj, = kn + rij,, k £ Z\{ 0}. (7) 


If the CFs behave as weakly interacting particles, they 
will form incompressible (topological) insulating states 
when filling an integer number of bands. Their band 
structure is given by a Harper-Hofstadter Hamiltonian 
with flux density n The densities n* at which filled 
bands are realized are therefore given by a Diophantine 
equation of the form ( 1 ) for the composite fermion sys¬ 
tem, n* = C*ri^ — D*, with integer parameters C*, D*. 
In composite fermion theory, one can explain the frac¬ 
tional QHE as an integer QHE of composite fermions 
[44], and for the lattice case one thus predicts incom¬ 
pressible states at n = n*. Hence, with (7), we find 


- D* 
kC* + 1 


> 0 . 


( 8 ) 


Choices of the parameters C* and D* for the composite 
fermion gap yield various candidates for incompressible 
states in the spectrum, given in terms of density, as il¬ 
lustrated previously as Fig. 1 in Ref. 14. 

In order to relate the densities ( 8 ) to the “filling fac¬ 
tor” of FQH systems there are several choices that can be 
made. One choice is to consider the ratio, = n/n^, 
of particle density to flux density, which is natural in 
the continuum limit, n$ —» 0 , where the bands of the 
Harper-Hofstadter model reduce to continuum Landau 
levels. More generally, and to allow connections to frac¬ 
tional Chern insulator models, the natural filling factor 
to consider is the ratio v = n/n s of the particle den¬ 
sity to the number density of single-particle states in the 
low-energy manifold (e.g. the lowest energy band if this 
manifold is a single band). We replace n$ in Eq. ( 8 ) via 
( 1 ) and obtain 


/l + fcC*\ D* n s D 
^ C* ) + C* ~ ~C + C' 


( 9 ) 
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In general, all parameters {k,C,C*, D, D*} contribute 
to determine eligible states. However, as we now de¬ 
scribe, for some important cases the ratios §? = § are 
equal, and the states can be characterized by a fixed 
filling factor v = n/n s . 

It is instructive to consider special cases. Firstly, the 
fractional quantum Hall states in a C = 1 band are recov¬ 
ered by choosing the manifold as the (lattice equivalent 
of the) lowest Landau level (C = 1, D = 0), and tak¬ 
ing general integers C* = t, D* = 0 for filling composite 
fermion states in the f-th Landau level. One recovers the 
usual Jain series of states v = t/(kt + 1). (In this case, 
n s = n 4>i so the two definitions of filling factor coincide.) 

Secondly, we can take both the CF bands and the ef¬ 
fective low-energy bands in the same subcell of the Hof- 
stadter spectrum, close to the flux = 1/|C|. As a 
concrete example, consider the gaps (5) and choose to 
fill r bands of composite fermions such that C* = rC , 
and D* = rsgn(C), and choose the lowest band with 
the given C, and D = sgn(C) for the manifold of single¬ 
particle states. We obtain states with filling factors [45] 


v 


C*—rC 


n 

n s 


r 

r\kC\ + 1 ’ 


r G Z\{0}, 


( 10 ) 


where states with r < 0 represent the generalization 
of negative flux attachment. In general, the low-energy 
manifold supporting these states will have many bands, 
but in the cases (4) this reduces to a single band. 

The sequence of filling factors (10), valid for any Chern 
number C ^ 0, is a core result of our Letter [46]. 

Several remarks are in order. The case with r — 1 can 
be seen as an analogue of the Laughlin state, in the sense 
that a single band of composite fermion states is filled. 
From the previous studies of the Laughlin state on the 
lattice in a C = 1 band, we can infer useful intuition 
on the likely stability of such states. The Laughlin state 
was shown to be stable up to flux densities ~ 0.4, 
i.e. it persists through 80% of the region in which a gap 
is open [47]. Likewise, the v = 2/3 hierarchy state was 
seen to be stable up to ~ 0.3 [14]. In the case of the 
states stabilized in subcells with \C\ > 1 bands, we note 
that the bands tend to have less dispersion, albeit maybe 
larger fluctuations in the band geometry. By analogy, the 
family of states (10) can be expected to be stable at a 
substantial distance from the respective cell boundaries. 

The reader will note that the prediction of compos¬ 
ite fermion theory (10) includes the Abelian states at 
filling factors v = l/(|fcC| + 1) that have recently been 
described in studies of Chern bands with Chern num¬ 
ber |C| > 1 [22-24], for both bosons and fermions, 
and which were described in terms of C flavour states 
[13, 25, 42, 43]. The derivation presented here demon¬ 
strates that these states are predicted also by the con¬ 
cept of flux attachment [11] leading to CF wave func¬ 
tions of the form described in Ref. [14]. While the C- 
flavour/multilayer language appears to require C copies 
of a C = 1 Brillouin zone, implying finite size geome¬ 
tries with a number of states N s mod C = 0, it was 
shown that a color-entangled formulation remedies this 
constraint [25]. Note that the hierarchy wave functions 


following from the CF construction [14] similarly do not 
require any constraint on the lattice geometry. 

The composite fermion theory makes a more general 
prediction, in that it does not require that the single par¬ 
ticle states making up the manifold n s are from a single 
band, as n$ can vary continuously in (9). Indeed, per¬ 
turbing a stable quantum liquid formed in a single-band 
configuration by an infinitesimal change in n$, the low- 
energy manifold splits up into (possibly infinitely) many 
bands, but we expect that the physics of the phase should 
be robust under this perturbation, providing a notion of 
adiabatic continuity that allows us to connect any band 
to the limit of the perfectly flat general Chern bands ob¬ 
tained as ri 0 —* 1/C [40], in line with the behaviour seen 
in C = 1 Harper-Hofstadter bands [48]. 

The filling factors (10) are analogous to the hierarchy 
states, which have been observed in Chern number one 
fractional Chern insulators [14, 49, 50]. Their proper¬ 
ties, in terms of quasiparticle charges and statistics were 
predicted by Kol and Read [11], as summarized in [40]. 
Unlike the lowest Landau level, \C\ > 1 bands support 
states with negative flux attachment (r < 0 ) even for 
|fc| = 1 , r = — 1 , so the corresponding series of states 
leads to novel filling factors. Numerical evidence for the 
v = 1 state in a C = 2 band obtained for Ckr = — 2 was 
provided by the current authors in Ref. [14], which is a 
special case in that it realizes an integer quantum Hall 
effect of bosons [51]. 

The limit of filling many CF Landau levels, r —> oo, 
which represents the equivalent of the half-filled Landau 
level, converges to 


lim v C '= rC 

r—t oo 


l 

WY 


( 11 ) 


At these points, the composite fermion spectrum resem¬ 
bles a Fermi-sea, as the band-gaps between the com¬ 
posite fermion levels decrease as 1/r and evolve into a 
quasi-continuum. In analogy to the half-filled continuum 
Landau levels, one may ask whether this filling can be 
susceptible to the equivalent of a CF pairing instabil¬ 
ity, or possibly more exotic states. In the C = 1 case, 
the possibility of a Moore-Read state at v = 1 is well 
known [52]. For the C = 2 band near = 1/2, a paired 
phase has been described in a related continuum model 
[43, 53], though the model does not provide a quantita¬ 
tive description of C = 2 bands [54]. If realized, it is 
expected that a paired phase at the fillings (11) would 
be non-Abelian in Chern bands with \C\ odd , while Ma- 
jorana quasiparticles will likely pair up for even Chern 
bands and thus recombine to yield an Abelian phase [53]. 

A case that has not yet been explored is the Abelian 
series of states (10). We examine the evidence for the 
presence of these composite fermion or hierarchy states 
on the basis of the band-projected Hamiltonian within 
the low-energy manifold, focusing on the single-band 
cases (4). The corresponding Hamiltonian, 77 pro j = 
T’lem’HT’lem, can be studied in the same framework as 
other fractional Chern insulator models [17]. The resid¬ 
ual dispersion of bands in the low-energy manifold could 
be of interest for studying phase transitions between frac- 
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FIG. 1. Spectra for Chern insulators of the series (10) in 
Harper-Hofstadter bands at flux densities from (4). For 
Chern number |C| = 2 we choose flux density n</, = 7/15, 
showing (a) v = 2/5 for IV = 10 (b) v = 2/3 for IV = 12 
bosons. For C = 3 we take = 7/20, and show (c) v = 2/7 
for TV = 8, (d) v = 2/5 for TV = 10. For lattice geometries, 
see legend. Multiplet structure of ground states shown in 
parentheses. 


tional quantum Hall liquids and condensed phases of 
bosons or Fermi liquid like states of fermions, respec¬ 
tively. However, here we choose to neglect the residual 
band dispersion, particularly as it vanishes quickly as 
» 1/|(7| [40]. Furthermore, we focus on the case of 
bosons with contact interactions V/, = USij. 

Our numerical study shows evidence supporting the 
existence of gapped quantum liquids at several filling 
factors of the series (10). Firstly, states are found for 
the cases r = 1 , \k\ = 1 , where the predictions of 
the filling factor v = l/(|Cj + 1) coincide between the 
composite fermion theory and the analyses in terms of 
Halperin multi-component [13, 42, 43] or color-entangled 
states [25] . The integer bosonic quantum Hall state with 
r = — 1 was discussed in Ref. [14]. Here, we present ev¬ 
idence for additional states, such as those with |r| = 2, 
with two filled composite fermion bands. In figure 1, we 
show their spectra for Chern bands with \C\ = 2 [pan¬ 
els a), b)], and C = |3| [panels c), d)]. All cases show 
the correct ground state degeneracies predicted by CF 
theory (g?qs = |1 + kC* |, see Refs. [11, 40]). The states 
with positive flux attachment, v = 2/5 (1(71 = 2) and 
tz = 2/7(|Cj = 3) have the clearest signature in terms of 
the magnitude of the gap to the average state spacing of 
excitations. The states with negative flux attachment, 
v = 2/3 (|<7| = 2) and v = 2/5 (|(7| = 3) also show a 
distinct separation of energy scales. 

A finite size scaling of the gap gives us further indica¬ 
tions of the stability of these phases. Figure 2 shows the 
gap scaling for several filling factors in the same \C\ = 2 
and \C\ = 3 bands as above. All systems we examined 
show the expected ground state degeneracy. In both 
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FIG. 2. Finite size scaling of the gap for states of the series 
(10), with r = +1, +2, —2. Left: a C = —2 band ( n$ = 7/15) 
for states at v = 1/3, v = 2/5, and v = 2/3. Right: Finite 
size scaling for (7 = 3 states (n^ = 7/20) at v = 1/4, v = 2/7, 
and v = 2/5. 


bands, the largest gap is found for the r = 1 state of 
the series (10), while the r = ±2 states have a slightly 
smaller gap in the finite size systems. The extrapolation 
to the thermodynamic limit is consistent with a non-zero 
gap for both the r = 1 and r = ±2 states for |(7| = 2. 
Data for the (7 = 3 cases are both noisier and include 
fewer system sizes. Nonetheless, the results are consis¬ 
tent with a non-zero gap. 

In the supplementary material, we briefly discuss par¬ 
ticle entanglement spectra of our states, as well as spec¬ 
tral flow under flux insertion [40]. In addition to these 
ground state properties, we have examined spectra un¬ 
der addition of “flux,” i.e., under changes of system size 
at fixed TV. We find low-lying bands of states consis¬ 
tent with an interpretation as quasiparticle states of an 
underlying quantum Hall liquid. We leave the detailed 
analysis of these features for a future publication. 

In conclusion, we have translated the composite- 
fermion or Chern-Simons theory into the language of 
fractional Chern insulators, leading to the prediction of 
a series of states with filling factors v = r/(r\kC\ + 1), 
for bosons (|fc| = 1) or fermions (|fc| = 2). This in¬ 
cludes, and provides an alternative description for, the 
series of states v = l/(|fc(7| + 1) that were observed 
in the literature on FCI for |(7| > 1. We have identi¬ 
fied flux densities where a single isolated band of Chern 
number (7 occurs at the bottom of the Hofstadter spec¬ 
trum. Finally, we have studied the many-body states 
of bosons with contact interactions under the projection 
into these Chern bands, identifying gapped states with 
the ground state degeneracies predicted by theory. While 
previous evidence had been given for the bosonic integer 
Chern insulator state with v = 1, r = —1, |(7| = 2 
[14], which was obtained for a hard-core interaction and 
without applying a band projection, the current results 
provide evidence for the wider applicability of compos¬ 
ite fermion theory, and its validity also for the band- 
projected Hamiltonian. Our results can be extended to 
general |(7| > 1 bands in other tight-binding models, 
and to the effective continuum limit n$ —» 1/1(7| via 
principles of adiabatic continuity. Further investigations 
should focus on the stability of fermionic states, the role 
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of long-range interactions and the detailed analysis of the 
ground states and excitations in terms of the composite 
fermion trial wave functions. 
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Supplemental Material 


A. SINGLE PARTICLE 
HARPER-HOFSTADTER HAMILTONIAN 

As stated in the main text, the Harper-Hofstadter 
Hamiltonian can be written in a gauge with explicit peri¬ 
odicity given by multiples of the lattice vectors {vi,v 2 } 
of an underlying periodic potential. To this end, we re¬ 
quire a pair of magnetic translation operators Tm [1] 
that respect the translations of the lattice, i.e., we seek 

ti = Tm(1 iVi), t 2 = Tm(^ 2 V 2 ), (A.l) 

and we choose U such that the area swept by the trans¬ 
lations encloses the smallest possible integer number of 
flux quanta N^. The cell can take geometries satisfying 

hhWi A v 2 | = Nffi, (A.2) 

where lil 2 = q and N£ = p yield the target flux density 
of = p/q. Solutions to (A.2) provide a “magnetic unit 
cell” (MUC) for the Hofstadter problem. 

The phases <pij in the single-particle Hamiltonian (3) 
[2] for a given gauge take the explicit form 

e f r ° 

4>ij = J A -dl + S^ij, (A.3) 

where 5(f>ij encodes the additional phase factor gener¬ 
ated by the magnetic translation operator (relative to 
canonical translations) for any hopping terms crossing 
the boundary of a MUC. For a given MUC with transla¬ 
tions {ijVj}, the entries tij exp(i(f>ij) of the single-particle 
tight-binding Hamiltonian are periodic under canoni¬ 
cal translations of the MUC. Hence, the many-body 
Hamiltonian can be simulated on finite-size tori of size 
N ce n = Lj x L 2 MUCs. For the Hofstadter model in the 
square geometry, we also denote indices i for , Li by the 
corresponding coordinate directions x, y. 


For our studies of the many-body Hamiltonian under 
band projection, we focus on the single-band cases (4) 
and retain N s = N ce n = L x x L y momentum sectors. 
Note that the total number of sites N s - lte = N x x N y = 
L x l x x Lyly N s , allowing for a large reduction of the 
corresponding Hilbert-space. The total particle density 
n s is small, hence band projection is expected to be valid 
even for strong repulsive interactions. Indeed, the states 
of the r = — 1 series were found to be stable for hard-core 
interactions in Ref. [3]. 


B. PROPERTIES OF HARPER-HOFSTADTER 
BANDS 

The Harper-Hofstadter bands are generally weakly dis¬ 
persive, particularly at the extremities of the spectrum. 
We display an example spectrum in Fig. S-3, for the case 
of 7i0 = 4/9 given in Harper’s gauge of an l x x l y = 9 x 1 
magnetic unit cell. On the scale of the entire spectrum, 
the dispersion is negligible (a), while detailed exami¬ 
nation reveals a residual dispersion, which carries a q- 
fold symmetry along the extended direction of the cho¬ 
sen magnetic unit cell (b). Owing to this high sym¬ 
metry, the Hofstadter spectrum is well suited to illus¬ 
trate a subtlety in calculating the Berry curvatures [4], 
two versions of which are shown in panels c) and d) for 
the lowest band. In the former case, we take Fourier- 
transforms with respect to the position of the unit cell, 
defining Bloch functions in reciprocal space via u a (k) = 
N~ 1/2 Er e lkR u(R + p a ), with sublattices labelled by a 
at sublattice offsets p a and origins of MUCs denoted by 
R. From these Bloch states, we obtain the first of the 
Berry curvatures, B\jq 1 shown in panel c), which breaks 
the (translational) symmetry of the problem. Taking the 
canonical Fourier transform with respect to the site po¬ 
sition, with u Q (k) = Nc l ^ 2 Er e' k ( R+p “)-u(R + p a ), we 
obtain the Berry curvature B shown in panel d), which 
manifestly respects the translational symmetry of the 
problem. These findings illustrate the general view that 
the position of sublattices in the unit cell must be taken 
into account in calculating the Berry curvature, as phys¬ 
ical properties probing B derive from the canonical po¬ 
sition operator [4]. 

In Fig. S-4, we demonstrate how the Berry curvature 
of the Harper-Hofstadter bands becomes perfectly flat 
near the points -A l/|Cj. We show that the extremal 
values of the Berry curvature approach the expected av¬ 
erage value exponentially in the parameter p of Eq. (5) 
in the main text. A similar exponential behavior is seen 
in the band dispersion. See also Ref. [5] for a discussion. 


C. ADIABATIC CONTINUITY OF QUANTUM 
LIQUIDS IN THE HOFSTADTER SPECTRUM 

There are two different kinds of adiabatic continu¬ 
ity which are relevant in the context of the Harper- 
Hofstadter model. 
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FIG. S-3. Single-particle properties for the Hofstadter Hamiltonian at flux-density = 4/9 for a l x X l y = 9 x 1 magnetic 
unit cell (q = 9): a) Dispersion for all nine bands, b) Detailed view for the dispersion of the lowest band, c) Berry curvature 
of the lowest band Hue, when Fourier transforms are taken relative to the unit cell (incorrect, see main text), d) canonical 
Berry curvature of the lowest band B. 


Firstly, we note that the composite fermion theory in¬ 
corporates an inherent notion of continuity for both n 
and ri 0 in the definition of preferred densities, and thus 
of the many-body spectrum under changes of these pa¬ 
rameters satisfying (9). The prediction is that a sys¬ 
tem in which the low-energy manifold consists of a sin¬ 
gle band should manifest the same topological order as 
a system nearby in parameter space, in which n$ is an 
irrational number, and this manifold fragments into in¬ 
finitely many bands. In this sense, band filling is not 
the most useful way to analyze the problem. Rather, 
the useful number of states n s to consider is contained 
within the (ensemble of) band(s) below a distinct gap 
in the Hofstadter spectrum. In particular, this sense of 
adiabatic continuity allows us to change the flux density 
gradually until A— 1 / 1 G7| is small, and the ef¬ 
fective magnetic length becomes large. This yields an 
effective continuum limit, providing the equivalent of a 
Landau-level with Chern number C, in which both the 
band dispersion and the band geometry are perfectly flat. 

Secondly, we can use the predictions of composite 


fermion theory to infer the behavior of general single¬ 
particle bands with Chern numbers \C\ > 1. In most 
models of Chern bands, we fix the number of sites per 
unit cell, and hence the number of bands as the ele¬ 
mentary input for the problem. Focusing on the cases 
where the low-energy manifold consists of a single band, 
the theory for the Hofstadter model predicts FQH states 
at the band fillings (10) in bands of Chern number Ch 
for the Hofstadter spectrum. As in the case of Chern 
number one bands [ 6 - 8 ], one can construct a sequence 
of many-body Hamiltonians that adiabatically deform 
the single-particle Hamiltonian of a suitable Hofstadter 
model to the corresponding Chern insulator model, as 
long as the Chern numbers for the relevant (occupied) 
bands are matched C = Ch- 

For the case of C = 1, we see that the adiabatic con¬ 
nection to the problem of continuum Landau levels [ 6 , 8 ] 
can be constructed also in the Hofstadter model, as the 
deformation within the gap C = 1, D = 0 yields the 
continuum limit for -4 0. The continuum limit is 
attractive as it provides a flat Berry curvature, which is 
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FIG. S-4. Deviation of the canonical Berry curvature from its 
average value, defined as A B = \B — Cq/(2n)\ at the maxima 
and minima of B within the Brillouin zone (in units of q/2n). 
Data is given for the lowest bands the flux densities according 
to the series of Eq. (4), with = p/q = p/{\C\p — sgnC). 
Note the logarithmic scale, so the decay is exponential as 
n^l/\C\. 


believed to be optimal for the stability of Chern insula¬ 
tors (see e.g., [4, 9, 10]). More generally, we can take 
effective continuum limits to perfectly flat |(7| > 1 bands 
at the points —> 1/|(7|. These points are analogous to 
the ricf, —> 0 case in terms of their large MUC, and their 
perfectly flat dispersion and Berry curvature. Note that 
a perfectly flat Berry curvature can also be obtained for 
finite size systems of square geometry [11], 

In order to construct the adiabatic continuation of a 
\C\ > 1 Chern bands of a generic tight-binding model to 
a Hofstadter model at finite A n^, there are two proven 
strategies. Firstly, one may establish a mapping between 
these models at the level of the single-particle orbitals, 
either in the Wannier basis [6, 12, 13] or more favourably 
in the Bloch basis [8] . Secondly, one may pursue a strat¬ 
egy of embedding the target model as a subset of a Hof¬ 
stadter lattice and then adiabatically tuning the respec¬ 
tive bonds [7]. 

D. GROUND-STATE DEGENERACY AND 
QUASIPARTICLE PROPERTIES 

For the sake of clarity, let us briefly re-state the re¬ 
sults by Kol and Read in our notation. Following their 
work [14], the Hall conductance a xy and quasiparticle 
charge e* h follow simply from considering a Laughlin- 
type thought experiment of flux insertion for the com¬ 
posite fermions. The Hall conductivity C* of composite 
fermions (the er™J of Kol and Read) follows from the ap¬ 
plication of the Streda formula to the composite fermion 
system. To transport a composite fermion, one needs to 
insert not just one flux quantum, but one also needs to 
supply kC* additional flux to compensate for the gauge- 
flux attached to the particles in (7). Again, from Laugh- 


lin’s argument the quasiparticle charge also follows, lead¬ 
ing to Kol and Read’s results [14] 


Gxy 


C* e 2 * e 

l + kC*~h’ gqh - 1 + kC* ‘ 


(D.l) 


Finally, their Chern-Simons theory yields the ground- 
state degeneracy g?gs and statistics angle d q h f° r the 
(Abelian) statistics of the quasiparticles as 


dGS 


|1 + kC*\, 


(iqh = 7T 



1 + kC*J 


(D.2) 


As an immediate consequence of these relations, we see 
that the many-body Chern number of the dcs-fold de¬ 
generate ground-state manifold is equal to the Chern 
number of the composite fermion gap C *. This matches 
the result obtained for a state with 1(7*1 = 2 in our pre¬ 
vious work [3]. 


E. PARTICLE-ENTANGLEMENT SPECTRA 

For the r = 1 candidates of our series (10), the particle 
entanglement spectra yield the known entanglement sec¬ 
tor count of states of the color-entangled Halperin states 
[15]. Note also that the entanglement spectra for (7 = 1 
bands of the Harper-Hofstadter model have been studied 
in Ref. [16]. 

The analysis of particle entanglement spectra for ad¬ 
ditional states with r / 1 for the series of hierarchy / 
composite fermion states (10) is difficult, as no simple 
counting rule for these spectra is known even for Chern 
number (7 = 1 bands. Correspondingly, these spectra 
have been little studied, and a detailed analysis of the 
particle entanglement spectra would be a complex topic 
on its own. We will not enter into details in the current 
paper, however, it is interesting to compare the structure 
of the PES at the same filling fractions in different Chern 
bands for the sake of illustration. In Fig. S-5, we con¬ 
sider the case of the entanglement spectrum for Na = 3 
out of N = 10 bosons at filling factor v = 2/3 as realized 
in the (7 = 2 Hofstadter bands at flux density = 7/15 
(a) in comparison to the PES in the (7 = 1 of the Ruby 
lattice model (b). The latter was previously discussed 
and shown in Fig. 6 of Ref. [17]. 

We find both similarities and differences between the 
corresponding spectrum of bosons at v = 2/3 in a (7 = 2 
band, Fig. S-5a), and the (7 = 1 data in Fig. S-5b). 
The overall scale of the entanglement energies is smaller 
for (7 = 2: we find the largest entanglement energies 
are £max — 10.7, compared to £max — 15.5 for (7=1. 
Correspondingly, the entanglement gap is also smaller, 
at A]? =2 ~ 1.8, compared to A]/ =1 ~ 4.5. However, the 
fraction of entanglement gap to the overall entanglement 
bandwidth S { = A^/(^ max — £ m ; n ) are quite similar with 
8^ =2 ~ 0.31, compared to ~ 0.42. 

For the (7=1 case, it was possible to infer the ex¬ 
pected number of states from the corresponding contin¬ 
uum FQHE in the lowest Landau level by virtue of a 
FCI to FQHE mapping [17]. By virtue of this property, 
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FIG. S-5. Particle entanglement spectra for composite fermion / hierarchy states of bosons at v = 2/3. Both panels show 
spectra at IV = 10 bosons, with Na = 3. We display the entanglement energy £ as a function of momentum sector K = 
k x L y + k y . a) PES for for the C = 2 band of the Harper-Hofstadter Hamiltonian at flux density = 7/15. b) Reference 
case of the PES in the C = 1 band of the Ruby lattice model, at the parameters chosen by Liu et al. [17]. Note the different 
scales on the respective j/-axes. Numbers in parentheses indicate the cumulative count of eigenstates from the bottom of the 
spectrum. The green dotted line indicates the principal entanglement gap. 


that work found that the count of entanglement eigen¬ 
values below the entanglement gap [indicated by a green 
line in our Fig. S-5b)[ is identical to the count inferred 
from the vanishing properties of the v = 2/3 compos¬ 
ite fermion state. In the case of higher Chern numbers, 
no corresponding continuum models are known, so we 
cannot take this correspondence as a given. For the ex¬ 
ample under consideration, we find that the number of 
entanglement states (as indicated in Fig. S-5) below the 
principal gap is the same in both the (7 = 2 and the 
(7=1 models, with 43 of a total 46 states in sectors 
with k x = 0 and 42 of 45 states in sectors with k x 0. 
However, the blocking of states below that gap is clearly 
different, so the lower-lying (7=1 entanglement gap 
identified with the Gaffnian state in [17] does not carry 
over. Other differences appear in the dependency of the 
entanglement count on the number of particles Na that 
are traced out from the density matrix. For example, 
Na = 4, we find 200 states below the largest entangle¬ 
ment gap for the (7=1 Ruby lattice model, while there 
are only 198 states for the (7 = 2 Hofstadter model. 

In conclusion, we find that the entanglement spectrum 
of the higher Chern number case bears a considerable 
resemblance to its (7 = 1 counterpart. Hence, the data 
suggests that the physical nature of the ground state is 
also the same in both cases, i.e. that of an incompressible 
quantum liquid. A detailed analysis would require the 
development of theoretical predictions for the entangle¬ 
ment sector count of the hierarchy wave functions, which 
could be based on, for example, the thin-torus framework 
of Ref. [18]. 


F. SPECTRAL FLOW 

The spectral response to the insertion of flux pro¬ 
vides a characteristic feature of incompressible quan¬ 
tum liquids. In our numerical simulations, we can 
probe this response by applying twisted boundary condi¬ 
tions to our periodic simulation cell, notably by setting 
boundary conditions for single-particle particle states as 
i />(r + L a e a ) = For a state in the fractional 

quantum Hall regime, we expect to see spectral flow 
among the states of the ground-state manifold, but no 
mixing among the latter and any higher excited states. 
Indeed, this scenario is realized in the candidates for frac¬ 
tional quantum Hall or FCI states at filling factors of the 
series (10). 

For illustration, we examine the spectral flow for states 
with |r| = 2 in a \C\ = 2 band in Fig. S-6, displaying 
filling factors zv = 2/5 and v = 2/3 for positive and 
negative flux attachment, respectively. We consider the 
same geometries shown in Fig. la,b) of the main text. In 
both these cases, all dcs ground states are found in the 
( k x , k v ) = (0, 0) sector, so these states are not protected 
by symmetry from hybridizing with each other. 

For the v = 2/5 state, we choose a system of N = 10 
particles N s ; te = 15 x 25 sites, in a gauge with magnetic 
unit cell of l x = 3 times l y = 5 sites (and N s = 25 states 
in the lowest band). We find that the (near-) degenerate 
ground-states cycle among each other, remaining well 
separated in energy from all excited states under varia¬ 
tion of (f> y . Considering near-degenerate points of pairs of 
ground-states as crossings, the ensemble of states returns 
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FIG. S-6. Spectral flow for FCI states in a \C\ =2 Hofstadter 
band at = 7/15. Note the splitting of the ground states is 
orders of magnitude smaller than the gap to the first excited 
state, so we choose to show the flow within the ground-state 
manifold only. For reference, the full spectrum at (j> x = <j> y = 
0 is shown to the left of the respective panels. Top: Spectral 
flow for N = 10 bosons in a C = —2 band at filling v = 2/5. 
Bottom: Spectral flow for N = 12 bosons in the same C = —2 
band at filling v = 2/3. Geometries are further described in 
the main text. 


to its initial ordering and energies only after an insertion 
of five units of flux, owing to the five-fold degeneracy. 

For the state at v = 2/3, we choose a system of N = 12 
particles in A7 site = 18 x 15 sites (N s = 18), using the 
same magnetic unit cell of 3 x 5 sites. A similar picture 


ensues, demonstrating a flow among the ground state 
energies with a 3 x 27r periodicity of the spectrum in cf> x . 

Like the other features that we have examined, these 
results support our identification of the partially filled 
Hofstadter bands at the fillings (10) as well-formed topo¬ 
logical quantum liquids. 
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